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Abstract. The strong CP problem was solved by Peccei & Quinn by introducing axions, which 
are a viable candidate for Dark Matter (DM). Here the PQ approach is modified so to yield also 
D ■ Dark Energy (DE), which arises in fair proportions, without tuning any extra parameter. DM and 

m ' DE arise from a single scalar field and, in the present ecpoch, are weakly coupled. Fluctuations have 

a fair evolution. The model is also fitted to WMAP first-year release, using a Markov Chain Monte 
■ Carlo technique, and performs as well as ACDM, coupled or uncoupled DE. Best-fit cosmological 

parameters for different models are mostly within 2-a level. The main peculiarity of the model is 
to favor high values of the Hubble parameter. 
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Most cosmological data, including Cosmic Microwave Background (CMB) 
anisotropics, Large Scale Structure (LSS), as well as data on SNIa [1], are fitted 
by a ACDM model with density parameters — 0.7, ~ 0.3, Q.}, ~ 0.04 (for 



2 ■ DE, non-relativistic matter and baryons), Hubble parameter h ~ 0.7 (in units of 100 

to '. km/s/Mpc) and primeval spectral index — 1. However, the success of such ACDM 

' does not hide its uneasiness. The parameters of standard CDM are still to be increased 

by one, in order to tune DE. Furthermore, if DE is ascribed to vacuum, this turns out 
^ , to be quite a fine tuning. This conceptual problem was eased by dynamical DE models 

^ I [2, 3]. They postulate the existence of an ad-hoc scalar field, self-interacting through a 

suitable effective potential, which depends, at least, on a further parameter (for SUGRA 
potential [3] here considered, this is an energy scale A or an exponent a). 

Within the frame of dynamical DE models, we tried to take a step forward [4]. Instead 
of invoking an ad-hoc interaction, we refer to the field introduced by Peccei & Quinn 
(PQ) [5] to solve the strong-CP problem. Such scheme was already shown to yield DM 
[6]. In [4] we slightly modify the PQ scheme, replacing the Nambu-Goldstone (NG) 
potential introduced ad-hoc, by a tracker potential [3] so to yield also DE. This scheme 
solves the strong-CP problem even more efficiently than the original PQ model. We 
shall call this cosmology dual-axion model. 

Its main peculiarity is naturally predicting DM-DE coupling. A number of authors 
discussed coupled DE models [7], where a parameter /3 fixes the coupling strength. 
Limits on /3 arise from linear analysis, comparing predictions with CMB [8] or, more 
efficiently, from non-linear analysis. Non linearity boosts the effects of coupling [9, 10]. 
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Dual-axion model has several advantages both in respect to ACDM and ordinary or 
coupled dynamical DE: (i) it requires no fine tuning; (ii) it adds no parameter to the 
standard PQ scheme. DM and DE arise in fair proportions by fixing the energy scale 
A of the tracker potential. Further, the model has no extra coupling parameter being 
the strength of the coupling set by the theory; (iii) it introduces no field or interaction, 
besides those required by particle physics. This scheme, however, leads to predictions 
(slightly) different from ACDM, for a number of observables. In principle, therefore, it 
can be falsified by data. 

The only degree of freedom still allowed is the choice of the tracker potential. Up 
to now, the dual-axion scheme has been explored just in association with a SUGRA 
potential. In this case, it predicts a fair growth of density fluctuations, so granting a 
viable picture for the LSS. 

The WMAP data on CMB anisotropics allow to submit the dual-axion model to 
further stringent tests, by comparing it with other cosmologies as ACDM, standard and 
interacting dynamical DE. None of the above models performs neatly better than the 
others. Apparently, the best fit is obtained by uncoupled DE with SUGRA potential. 



A SINGLE SCALAR FIELD TO ACCOUNT FOR DM AND DE 

The solutions of the strong CP problem proposed by PQ leads to one of the accepted 
models of DM. PQ suggested that d parameter, in the effective lagrangian term 

^e = ^eG-G (1) 

{as', strong coupling constant, G: gluon field tensor), causing CP violations in strong 
interactions, is a dynamical variable. Under suitable conditions d approaches zero in 
our epoch, so that the term (1) is suppressed, while residual Q oscillations yield DM [6]. 

By adding to the Standard Model a global U(l) symmetry which is spontaneously 
broken at a scale FpQ, 9 is then the phase of a complex field <i> = (j)e^^/V2 which, 
falling into one of the degenerated minima of an NG potential 

vm)=m^\'-F^Q?, (2) 

develops a vacuum expectation value < >= Fpg. The CP- violating term, arising 
around quark-hadron transition when qq condensates break the chiral symmetry, reads 



Vi = 



Y,{OiT)\qq\0{T))m, 
1 



(1-COS0) (3) 



(L^ extends over all quarks), so that 6 is no longer arbitrary, but shall be ruled by a 
suitable equation of motion. The term in square brackets, at T ~ 0, approaches m|/| 
(nij^ and fn'- TT-meson mass and decay constant). 

In the next Sections, we will discuss the work in [4], where the NG potential (2) is 
replaced by a tracker potential. Then, instead of settling on a value Fpq, continues to 
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evolve over cosmological times, at any T. As in the PQ case, the potential involves a 
complex field 4> and isU{l) invariant. 

At variance from the PQ case, however, the 6 evolution starts and continues while 
also ^ is still evolving. This goes on until our epoch, when <j> is expected to account 
for DE, while, superimposed to such slow evolution, faster transversal 6 oscillations 
occur, accounting for DM. As it can be expected, however, DM and DE are dynamically 
coupled, although this coupling weakens as we approach the present era. 

Lagrangian theory 

In the dual-axion model we start from the lagrangian jSf = ^/^{gny^n^^v^ — 
y(|^»|)} which can be rewritten in terms of and 6, adding also the term breaking 
the symmetry, as follows: 

^ = ^/^{{l/2)g^,y[^^(|)^y(|) + (|)%^^ye]-V{(|>)-m\T,(j))(j)\l-cose)} . (4) 

Here guv is the metric tensor. We shall assume that ds^ = g^ydx^dx^ = a^{dx^ — 
r\ijdxidxj), so that a is the scale factor, t is the conformal time; greek (latin) indeces run 
from to 3 (1 to 3); dots indicate differentiation with respect to T. The mass behavior 
for T ~ AgcD will be detailed in the next Section. The equations of motion read 

+ 2(a/a + 0/0) 0+mVsin0 = , (5) 

<f> + 2{d/a)^+a^V'{^) = ^d^, (6) 

while the energy densities pe,^ = Pe,(i>;km + Pe,^-pot and pressures pe^^ = pQ,^;kin - 
Pe,(ti;pot-, under the condition ^ 1, are obtainable from 

Pe,kin = (l>^d^/2a^ , Pe,pot = m\T, (^)(^^e^/2 , 

P<!),kin = (j)^/2a^ , p^^pot = V{(j)) . (7) 

Axion mass 

According to eq. (5), the axion field begins to oscillate when m{T,^)a ~ 2(d/a + 
0/0 ). In the dual-axion model, just as for PQ, axions become massive when the chiral 
symmetry is broken by the formation of the qq condensate at T ~ Aqcd- Around such 
T, therefore, the axion mass grows rapidly. In the dual-axion model, however, a slower 
growth takes place also later, because of the evolution of 0. Then m{T, 0) is 

mo(0) = m^f^/(l) = (0.0062/0) GeV . (8) 

Since ~ today, the axion mass is now mo ~ 5 • 10~^^eV, while, according to [11], 
at high T : 

m(r,0)~O.lmo(0)(A2cz)/r)3-^ . (9) 
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This expression must be interpolated with eq. (8), to study the fluctuation onset for 
T ~ Aqcd- Details on interpolation can be found in [4]. 



The case of SUGRA potential 

When performs many oscillations within a Hubble time, then (pe^kin) — {Pe,pot) 
and (pq) ~ 0. By using eqs. (5),(6),(7), it is easy to see that 

m , , m 

pe+3Hp0 = -pe , p^+3H{p^+p^) = — pe . (10) 
m m 

When mis given by Eq (9) ,m/m = — ^/<^ — 3.8 T/r. At T ~ 0, instead, m/m~ —^l^- 
Here below, the indices q, ^ will be replaced by dm,de- Eqs. (10) show an energy 
exchange between DM and DE. The former eq. (10) can then be integrated, yielding 
pDM °^m./a^. This law holds also when T <^ Aqcd, and then the usual behavior p£,M °^ 
is modified, becoming Pdm<^^^ — const. Let us now assume that the potential reads 

V(0) = (A«+V0«)exp(4;r0V'"2) (j^) 

(no 6 dependence); in the radiative era, it will then be 0"+^ = gaA^'^^a^z^ , with 
gee = a(a + 2)^/4(a + 6). This tracker solution holds until we approach the quark- 
hadron transition. Then, in Eq. (6), the DE-DM coupling term, ^6^, exceeds a^V' and 
we enter a different (tracking) regime. This is shown in detail in the left plot (panel 
(a)) of Fig.(l) obtained for D.m — 0.3, D-i, — 0.03 and Hubble constant h — 0.7 (in units 
of ICQ km/s/Mpc). Panel (b) then shows the low-z behavior, since DE energy density 
exceeds radiation (z — 100) and then overcomes baryons (z— 10) and DM (z ~ 3). In 
the central plot of Fig.(l) a landscape picture for all energy densities p,- (/ = r, b, 0, ^, 
i.e. radiation, baryons, DM, DE), down to a = 1 (top panel) and the related behaviors of 
the density parameters (bottom panel) are shown. 

In general, once Q.de (at z = 0) is assigned, a model with dynamical (coupled or 
uncoupled) DE is not yet univocally determined. For instance, the potential (11) depends 
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TABLE 1. SUGRA parameters for uncoupled DE (left), constant coupling (center) and (p model 
(right). 



uncoupled SUGRA 



Q.bh^ 0.025 0.001 





0.12 


0.02 


h 


0.63 


0.06 


T 


0.21 


0.07 


ris 


1.04 


0.04 


A 


0.97 


0.13 


X 


3.0 


7.7 



SUGRA with C=const 

X (x) 



Q.bh^ 0.024 0.001 





0.11 


0.02 


h 


0.74 


0.11 


T 


0.18 


0.07 


n.s- 


1.03 


0.04 


A 


0.92 


0.14 


X 


-0.5 


7.6 


P 


0.10 


0.07 



SUGRA with C=(;) i 

X (x) 



Q.bh^ 0.025 0.001 





0.11 


0.02 


h 


0.93 


0.05 


T 


0.26 


0.04 


Us 


1.23 


0.04 


A 


1.17 


0.10 


X 


4.8 


2.4 



on the parameters a and A and one of them can still be arbitrarily fixed. In dual-axion 
model such arbitrariness no longer exists. The observational value of the densities in 
the world forces the scale factor when oscillations start to lay about the quark- 
hadron transition, while also A is substantially fixed. For Q-dm = 0-27 we obtain A ~ 
1.5 ■ lOi^GeV and an ~ lO^^l But, when D.dm goes from 0.2 to 0.4, logio(A/GeV) 
(almost) linearly runs in the narrow interval 10.05-10.39, while a/, steadily lays at the 
eve of the quark-hadron transition. For more details on this point see [4]. 



Evolution of inhomogeneities 

Besides of predicting fair ratios between the world components, a viable model should 
also allow the formation of structures in the world. 

The dual-axion model belongs to the class of coupled DE models treated by Amen- 
dola [12], with a time dependent coupling C ( ) = 1 / . Fluctuations evolution is then ob- 
tained by solving the equations in [12], with the above C(0). The behavior shown in the 
right plot of Fig.( 1) is then found. The bottom panel compares DM fluctuation evolutions 
in the the dual-axion model (solid curves), with those in an analogous ACDM model 
(long-dashed curves) and in a coupled DE model with constant coupling C ~ {C{^)) 
(short-dashed curves). As shown by the plots, the overall growth, from recombination 
to now is similar in dual-axion and ACDM models, being quite smaller than in DE mod- 
els with constant coupling. The differences of dual-axion from ACDM are: (i) objects 
form earlier and (ii) baryon fluctuations keep below DM fluctuations until very recently. 



COIVIPARISON WITH WIVIAP DATA 

We have tested the dual-axion model against CMB data, together with other dynamical 
DE cosmologies. We used the Markov Chain Monte Carlo program [13] also used in the 
original analysis of WMAP first-year data [14] to constrain a flat ACDM in respect to 
six parameters: Clhh^, Q-mh^, K ris, the fluctuation amplitude A and the optical depth x. 
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FIGURE 2. Cj spectra for the best fit SUGRA (solid line), constant coupling (dotted line), 0"'- 
coupling (dashed) and dual-axion (dot-dashed) models. 



In our analysis, three classes of DE were considered: (i) uncoupled SUGRA DE, 
requiring an extra parameter A = logi()(A/GeV), the energy scale in the potential (12). 
(ii) Constant coupling SUGRA DE, with a further parameter /3 °« C. (iii) Coupled models 
with C = keeping just A as a free parameter. The (iii) class of model includes the 
dual-axion case, for which, however A is set by the requirement that Q.oe lays in a 
fair range so that A and Q.de, are no longer independent. Then, we tested whether data 
constrain A into a fair region, turning a generic model into a dual-axion model. 

The basic results are summarized in the Table 1 . For each model we list the expecta- 
tion values {x) of each parameter x and the related variance o^- 

A first point worth outlining is that SUGRA (uncoupled) models, bearing precise 
advantages in respect to ACDM, are consistent with WMAP data. 

In uncoupled or costant coupling SUGRA models opacity (t) is pushed to values even 
greater than in ACDM [15]. Greater t's have an indirect impact also on Q.hh^ whose 
best-fit value becomes greater, although consistent with ACDM within 1-C7. 

Parameters are more strongly constrained in models. In particular, WMAP data 
yield constraints on A for models and the 2-a A-interval ranges from ~ 10 to 
~ 3 ■ lO^^GeV, so including the dual-axion model. The main puzzling feature of 
models is that large h is favored: the best-fit 2-a interval does not extend below 0.85 . 
The problem is more severe for dual-axion A 's values. This model tends to displace the 
first Cf peak to greater I (smaller angular scales) as coupling does in any case does. 
The model, however, has no extra coupling parameter and the intensity of coupling is 
controlled by the scale A. Increasing this scale requires a more effective compensation 
and greater values of h are favored. This effect appears related to the choice of SUGRA 
potentials, which is just meant to provide a concrete framework for the dual axion model. 

The fits to WMAP data yield similar for all models ranging from 1.064 (no 
coupling) to 1.074 (0^' coupling) [4]. Fig. (2) compares the Cf spectra for all best- 
fit models (apart of ACDM). At large / all models yield similar behaviors and this is 
why no model category prevails. 
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CONCLUSIONS 



Axions have been a good candidate for DM since the late Seventies. They arise from the 
solution of the strong CP problem proposed by PQ. Here we showed that their model 
has a simple and natural generalization that also yields DE adding no parameter to the 
standard PQ scheme. A complex scalar fields 4>, arising in the solution of the strong CP 
problem, accounts for both DE and DM: as in the PQ model, in eq. (1) d is turned into 
a dynamical variable, the phase of 4>. Here, however, instead of taking a constant value, 
|4>| increases in time, approaching mp by our cosmic epoch, when it is DE; meanwhile, 
is driven to approach zero, still performing harmonic oscillations which are axion 
DM. The critical scale factor when oscillations start, lays at the eve of the quark- 
hadron transition, because of the rapid increase of the axion mass m(r, 0). The scale a;, 
is essentially model independent and no appreciable displacement can be expected just 
varying Q.dm while rather high values of h tend to be preferred. This implies that the 
scale A in the SUGRA potential is almost model independent [4]. Therefore, the unique 
setting of Q-dm fixes A also providing DM and DE in fair proportions and simultaneously 
solving the strong CP problem. 

The fits to WMAP first-year data of uncoupled, constant-coupling and SUGRA 
models yield similar ;t^'s. At variance from other models, in models, CMB data 
constrain A and it is significant that the allowed range includes values consistent with 
the dual-axion model. 
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